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Abstract 

We discuss exact results for the full nonperturbative effective superpotentials of four 
dimensional A/" = 1 supersymmetric U{N) gauge theories with additional chiral su- 
perfield in the adjoint representation and the free energies of the related zero dimen- 
sional bosonic matrix models with polynomial potentials in the planar limit using 
the Dijkgraaf-Vafa matrix model prescription and integrating in and out. The exact 
effective superpotentials are produced including the leading Veneziano-Yankielowicz 
term directly from the matrix models. We also discuss how to use integrating in and 
out as a tool to do random matrix integrals in the large limit. 
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1 Introduction 



The exact nonperturbative effective superpotential of A/" = 2 supersymmetric gauge tlieories witfi all 
instanton corrections is completely determined witli a knowledge of the perturbative UV physics, the 
singularities in moduli space and the monodromies around the singularities, The field content of 
a pure Af = 2 gauge theory is a vector field A^, two Weyl fermions Aq, and ipa, and a complex scalar 
A. In A/" = 1 language, these are a combination of a field strength chiral superfield Wa containing 
Aq, and A^, and a scalar chiral superfield $ containing A and ipa all transforming in the adjoint 
representation. Suppose this pure M = 2 theory is perturbed by a tree level superpotential including 
a mass term for $. When one integrates out $, the low energy theory below the mass of $ reduces 
to A/" = 1. About a decade ago, various examples of exact superpotentials of AA = 1 supersymmetric 
gauge theories were obtained making use of holomorphy and symmetry arguments. See ^ for a 
review. However, for such general tree level perturbations, combinations of parameters which are 
not protected by symmetries appear. These parameters can come in the effective superpotential 
with any degree. Dijkgraaf and Vafa [2j found that the effective superpotentials of A/" = 1 theories 
obtained by such deformations oi M = 2 could be computed by using matrix models. 

Consider pure A/" = 1 supersymmetric U{N) gauge theory. The SU{N) subgroup confines in 
the infrared and using the one loop running of the gauge coupling coefficient the Lagrangian that 
describes the theory at the cut off can be written as 



f (f e 3N \og{^) S + C.C., (1.1) 
J Ale 



where 5* is the glueball superfield defined in terms of the gauge chiral superfield as 

S=--^TtW^Wo., (1.2) 

Ai is the scale of the A/" = 1 SU{N) gauge theory and Ai^ is the UV cutoff. The point is that 
from p.l|) one sees that 3A^log(Ai) linearly couples to S and acts as its source. On the other hand, 
gaugino condensation gives the nonperturbative effective superpotential 

W,s = NAl (1.3) 

Note that we have suppressed the phases e^'^^^^^ , k = to N — 1, in p.3|) due to the breaking 
of the Z2N R symmetry down to Z2 by gaugino condensation. We can integrate in 5* to p.3|) 
with 3A^log(Ai) as its source and calculate the glueball superpotential while still keeping Ai as a 
parameter in the theory by introducing an auxiliary field A and minimizing 

W = NA^ -3NS\og{^) (1.4) 

Ai 

with A which gives the Veneziano-Yankielowicz effective superpotential [7j 

W,s = NS-NS\og{^l (1.5) 
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Integrating out S in ()1.4p simply gives back the original nonperturbative superpotential p.3|) . 

Suppose one adds a chiral superfield $ in the adjoint representation with a tree level superpo- 
tential 

n 

mree = V^Tr$^+^ (1.6) 

^ p+l 

p=l 

to the Af = 1 U (N) gauge theory. For instance, for the specific case of a cubic tree level super- 
potential with n = 2 in ()1.6p and the gauge symmetry in the low energy theory unbroken, there 
is a combination of parameters glS/g2 that has no charge under all symmetries. This parameter 
can appear with any power and the instanton corrections to the effective nonperturbative superpo- 
tential can be written as YlV=i^k{.glS / g\Y S where are constants. Our interest is to find exact 
analytic expressions for the glueball effective superpotential whose series expansion produces all 
terms including the leading Veneziano-Yankielowicz and the instanton corrections and to compute 
the free energies of the corresponding zero dimensional bosonic matrix models using integrating in 
jHj and out techniques. The relation between the Dijkgraaf-Vafa matrix model prescription and 
the integrating in method for U [N) gauge theories was discussed in IHj and our approach on the 
integrating in side here will be along similar lines. 

For completeness we will start with a very brief review of a zero dimensional U{N) bosonic 
matrix model as originally developed in [Hj and reviewed in pj. We will then compute the complete 
exact nonperturbative superpotentials for a quadratic, cubic and quartic tree level superpotentials 
using the Dijkgraaf-Vafa matrix model prescription. The free energies in the matrix model for cubic 
and quartic potentials were originally computed in [Hj and the corresponding instanton corrections to 
effective superpotentials were computed in jTUllII] using the Dijkgraaf-Vafa prescription. However, 
the leading Veneziano-Yankielowicz term was added by hand. Here we will point out the proper 
normalization of the partition function and scheme that produces the exact effective superpotential 
including the Veneziano-Yankielowicz term directly from the matrix model. We will then continue 
with computing the full exact effective superpotential for quadratic, cubic and quartic tree level 
deformations using integrating in and out. The results from the two approaches exactly agree to 
all order including the leading Veneziano-Yankielowicz term and the instanton corrections. We will 
also provide a scheme for using integrating in and out to do random matrix integrals. For the case 
of one-cut U{N) bosonic matrix models we will discuss in this note, integrating in and out combined 
with the Dijkgraaf-Vafa matrix model prescription provides a very simple tool to do random matrix 
integrals in the planar limit for any tree level polynomial potential. 

2 Matrix model 

Consider a zero dimensional U{N) bosonic matrix model with a potential given by ()1.6p . $ is an 
N X N matrix. Extremizing the potential generically gives n distinct values a/ of $, where / = 1 to 
n. If $ is taken to have classical eigenvalues aj each with degeneracy Nj such that = ^^1=1 
then the gauge symmetry in the low energy theory is broken to n/=i U{Ni) . Our interest in this 
note is the case where all aj are equal and the gauge symmetry is preserved. We will take aj = 0. 
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Now consider the partition function 



c^$e"^^*^'=^ (2.1^ 



Vo\{U{N)) 

The large N (or small Qg) limit is done with N ^ oo and (^s — ^ where 

S = Ngs (2.2) 

is fixed. The perturbative expansion of (j2.ip leads to a sum of Feynman diagrams over Riemann 
surfaces. In the large N limit only the leading order sum over a genus zero planar surface survives 
and the partition function becomes 

Z = e~^, (2.3) 

where jFg is the free energy in the planar limit. 

In the large N limit, the integral is done using the standard matrix model technology where one 
starts with transforming the integral from $ to a set of eigenvalues Aj of $. The matrix integral 
(12.111 when transformed to integral over eigenvalues Aj becomes 

Z= f JJrfA.e-^S«^("')+^^-^^°^''-'^l, (2.4) 

i 

where we write the tree level superpotential in (jl.fjj) as Wtree{^) = Tr H^($). The indices i and j in 
this section run over 1 to A^. The second term in the exponent of the integrand in ()2.4|) comes from 
the Jacobian in the transformation of the measure of the integral from the independent components 
of $ to the eigenvalues. This term gives rise to a repulsive interaction between the eigenvalues 
resulting in distributions of eigenvalues around each extremum point. In the case where the original 
gauge symmetry is preserved in the low energy theory and the extremum point is taken to be at 
the origin, all the eigenvalues are distributed around the origin. The equation of motion of the 
eigenvalues is 

ivr'(AO -2^-^ = 0. (2.5) 
Equation (j2.5j) is difficult to solve directly in practice and one introduces the resolvent 

i 

where x is complex. Note that w[x) has the large x asymptotic behavior 

w{x) -l/x. (2.7) 

The equation of motion of the eigenvalues in the large limit gives a quadratic equation for the 
resolvent, 

w{xf + \w'{x)w{x) + ^f{x) = (2.8) 



S ' ' ' ' AS 
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where 

i 

is a polynomial function of degree n — 1. ()2.8p has the solution 

wix) = (W'ix) ± ^/W'ixy - /(x)) . (2.10) 

Because W^'(x) and /(x) are polynomial functions, any singular behavior of w{x) comes from the 
square root in the second term in ()2.1()j) . Note that W^'(x)^ — f{x) is a polynomial of degree 2n 
and thus w{x) has in general the same number of branch points. This singularity structure and the 
asymptotic large x behavior of w{x) are enough to completely determine w{x) in the one branch 
cut case where n — 1 roots of W'{xy — f{x) come in pair. In the large N limit, one defines the 
eigenvalue density, 

'^(^) = i^E^(^^-^)' (2.11) 

i 

such that J p{X)dX = 1. The point is that w{x) and p{x) are related, 

,2,12) 

Therefore, once w{x) is obtained, we can invert (j2.1H) to compute the eigenvalue density, 

p{X) = ^(w{\ + iO) -w{\-iO)). (2.13) 

With p(A) in hand, the free energy in the large N limit follows from ()2.3j) and ()2.4|1 . 

^o = S J d\p{\)W{\) j J d\dpp{\)p{p) log |A - (2.14) 

Using the equations of motion given by ()2.5p in the large N limit, ()2.14j) can be written in a more 
convenient form, 

J^o = ]^S j dXp{\)Wi\)- j dAp(A)log|A|. (2.15) 
Following Dijkgraaf-Vafa 2J, the effective superpotential is given by 

W,s = N^ + NSlog{A'), (2.16) 

where A^log(A^) is a source for S and A is the nonperturbative scale in the M = 2 gauge theory. 

The important point and scheme we want to point out is that with the proper normalization 
of the partition function as given by ()2.H) and the planar limit of the free energy that follows 
from ()2.4j) . the exact effective superpotential is produced using ()2.1fij) without a need for adding 
the Veneziano-Yankielowicz term by hand. We will explicitly see in subsection (12.11) that the 
exact Veneziano-Yankielowicz superpotential follows from a quadratic tree level potential with this 
scheme. We will also derive the appropriate free energies in the planar limit for cubic and quartic 
tree level superpotentials and show that the full exact effective superpotentials including both the 
Veneziano-Yankielowicz term and all instanton correction are produced. 
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2.1 Quadratic potential 

Consider a quadratic mass term with n = 1 and g2 = mm Wtree given by Now W{x) = 

has only one extremum point and it occurs at the origin. The most general form of w{x) follows 

from (jfll) and dTTIUl . 

w{x) = —^{mx — m\/ x^ — a^), (2.17) 
where a = 2^ Sjm. The Wigner semicircle density distribution of the eigenvalues follows from 

fiTTTji in dnsi), 

p(A) = (2.18) 

for A G [—a, a]. Note that because VFtree is even, the expectation value of $ vanishes and U{N^ 
automatically reduces to SU{N). Using ()2.18|) and W{X) = mA^/2 in ()2.15|) . we obtain the free 
energy 

:Fo = ls^-lsHog{^). (2.19) 
4 2 m 

The effective superpotential follows from ()2.19p in ()2.16p . 

W,s = NS-NSlogi^\ (2.20) 

where Ai = (A^m)^/^ is the scale of the lower energy pure Af = 1 theory related to the scale A of 
the higher energy Af = 2 theory by threshold matching of the gauge coupling running at the energy 
scale m. Now ()2.20|) is exactly the Veneziano-Yankielowicz superpotential. It was discussed in j^lT^ 
that the volume of the U{N) gauge group gives a contribution like ^S*^ log(S') to the free energy in 
the planar limit. For a recent discussion on this point and on adding the Veneziano-Yankielowicz 
term by hand, see JT] . What we have shown here is that with the proper normalization and scheme 
discussed earlier in this section, the appropriate free energy and the exact Veneziano-Yankielowicz 
superpotential are obtained. There is no need for adding the Veneziano-Yankielowicz term by hand 
in computing the exact effective superpotential for any tree level superpotential deformation using 
this scheme. 



2.2 Cubic potential 

Next consider the cubic tree level superpotential, 

Wtree = ImTr ($2) + Tr ($3). (2.21) 

Analytic expressions for the free energy of this model and the quartic potential in the next subsection 
were first obtained in |H] and the instanton corrections in the effective superpotential using the 
Dijkgraaf-Vafa prescription have been computed in |10^ with the Veneziano-Yankielowicz term added 
by hand. Here we will compute the effective superpotential including the Veneziano-Yankielowicz 
term directly. Now there are two critical points, one at a; = and the other at a; = —m/g. For the 
case where the U{N) gauge symmetry is unbroken, we need to consider only one critical point and 
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we choose the one at the origin. In the quantum theory, all eignevalues are then distributed in an 
interval [a, b] enclosing the extremum point at a; = 0. The resolvent is then completely determined 
noting that W'{xy — f{x) in ()2.1()j) has one double root and using the asymptotic behavior ()2.7|1 . 
and it is given by 



^(9^'^ + ma; - {gx + ^(a + b)g + m)^{x - a){x - b)^, (2.22) 



where 



-g{a + by + -m{a + b) + -g{b - af = (2.23) 

and 

{b- aY[{a + b)g + m\-l6S = 0. (2.24) 
Combining the two conditions ()2.23p and ()2.24|) and defining 

^ = ^(« + &), (2.25) 

we obtain 

2o2 

a{l+(T){l + 2a) + ^S = (2.26) 

with solution 
where 



^=(./432«^-l-12^«4)'". (2.28) 



We have chosen the solution a to the cubic equation ()2.2(jj) such that a and b are real and the 
eigenvalues are distributed on the real axis of x for real g'^S/m?. The eigenvalue density follows 
from ()2.22|) and ()2.13|) . and writing it in terms of a and making a change of variable, 

P^y^ = 2^ (^^ + + ^''^'") V^o"^'' ^2-29) 
where y = X — ma/g and yo = 2a/ S/ ((1 + 2a)m). The eigenvalues are distributed in the range 
A G [a, b] or equivalently in y E [—yo, yo]- 

The free energy can then be calculated using ()2.29|) in ()2.15|) which gives 

= -2^ ''^^WT^)^ ^ m + a)(l + 2aY " ^^'''''^ 
Using (ITT^ in (ITT!)|l and noting 



da _ o-(l + (t)(1 + 2a) 
as ~ 6(t2 + 6a + 1 



(2.31) 



from ()2.26j) . we obtain 



where again Ai = (A^m)^/^ is the scale of the low energy A/" = 1 theory. We will make a series 
expansion of (|2.32|) in subsection (|3.2p and see that it produces the exact effective superpotential 
including the Veneziano-Yankielowicz term and the instanton corrections to all order. 
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2.3 Quartic potential 

Finally, let us consider the quartic superpotential 

l^tree = ^mTr($2) + ^(^Tr (<D^). (2.33) 

In this case, again for the one-cut case in which the gauge symmetry is preserved and all eigenvalues 
are distributed around the origin, the resolvent is computed noting that H^'(x)^ — /(x) in ()2.1()j) 
has two double roots and using the asymptotic behavior ()2.7|) . 

w{x) = (gx^ + mx — {gx^ H — ga^ + ra)\/ x'^ — a?^ , (2.34) 

where 

?,ga^ + Ama^ - 165 = 0. (2.35) 
The eigenvalue density follows from (j2.34|) in ()2.13|) 



P(A) = ^{gX"" + -ga' + m)v/^^^, (2.36) 



Using ()2.33j) and (j2.36|l in ()2.15|1 . we obtain the free energy 



= -\s''\og{^) + ^i-m'a^ + A{)mSa^ + 144^2). (2.37) 



2 °M ' 384 

Noting from that 

da 



dS 3ga^ + 2ma ' 
the effective superpotential is obtained using ()2.37p in ()2.16|) . 



^ 4 ^ 72c/a4 + 48ma2 48 ' ^ ' 



where 



a^ = |^(Jl + 124-1). (2.40) 
Sg V wi^ 



(2.38) 



We will show in subsection ()3.3|) that making a series expansion of ()2.39p gives the correct leading 
order Veneziano-Yankielowicz term in addition to all instanton corrections. 



3 Integrating in and out 

The theory without the tree level deformation is simply pure Af = 2 U{N) gauge theory and it 
has been well studied. The quantum moduli space of this theory can be parameterized by the 
hyperelliptic curve jH US] 

l/2 = P^(x)2-4A2^ (3.1) 
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where y and x are complex coordinated such that 



AT 



PN{x) = det{xl-^)=W{x-Xi), (3.2) 



i=l 



Xi being the diagonal elements of $, which can be made diagonal using D- flatness conditions. A 
gauge invariant parametrization of the moduli space can be given by coordinates 



N 



u^ = Ti^ = ^Xi (3.3) 



i=l 

and 

N 



u„ = -Tr ($ + uJNY = - \"(xi + Ui/NY, p = 2toN. (3.4) 
p p 

The quantum moduli space given by p.ip describes a genus — 1 Riemann surface with two 
types of cycles. One type (/3-cycle) is related to the handles of the surface and the second type 
(a-cycle) connects different handles. A dual scalar field is given by integral of a meromorphic 
differential over an a-cycle. A dual scalar field vanishes and the associated monopole (or dyon) 
becomes massless when an a-cycle degenerates. The strong coupling singularities of the quantum 
moduli space correspond to regions where monopoles become massless. Our interest here is the 
case where the U{N) gauge symmetry is preserved. This corresponds to the case in which all the 
a-cycles vanish and the genus A^ — 1 Riemann surface degenerates into a sphere. The effective field 
theory which describes the low energy excitations near these singularities is described in terms of 
the dual scalar and monopole fields and the exact effective superpotential is given by, following the 
nonrenormalization and linearity principle of [H], 

N-l n+1 
W =^ AD,ni{u)EmEm + QpUp, (3.5) 

m=l p=2 

where and E^ denote the monopole chiral multiplet which becomes massless at the m**^ singu- 
larity and are associated with the dual scalar field A^^rn- The second term in ()3.5j) is the tree level 
deformation. The steps to compute the effective glueball superpotential are integrating out Em, 
Em, Up and Ui and integrating in S. The equations of motion obtained by minimizing ()3.5p with 

Em, Em and Up are 

AdAu) = 0, (3.6) 



N~l 



2^ ^m^m + 5-^ = 0. (3.7) 



m=l 



Putting ()3.(ij] in ()3.5j) gives 



n+1 

W^s = {^gpUp)\AoM^)=o- (3.8) 

p=2 
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Now ()3.7|) shows the standard confinement of monopoles and ()3.6|) expresses the vanishing of the dual 
scalar fields and the masslessness of the monopoles at the singularities. The effective superpotential 
is thus given by the tree level deformation evaluated with the constraint that all the dual scalar 
fields vanish. The glueball effective superpotential is then obtained by integrating out Ui and 
integrating in S while at the same time keeping the nonperturbative scale as a parameter in the 
leading Veneziano-Yankielowicz part of the superpotential. This can be done by introducing an 
auxiliary field A and minimizing 

^ A 

J2(9pUp)\AoMu)=o (with A ^A)- 2iV51og(-). (3.9) 

p=2 

with A and Ui. Note that when all the a-cycles simultaneously degenerate, P^(a;)2-4A2^ has two 
single roots and — 1 double roots and the hyperelliptic curve factorizes. ^1 IH] -Pat is then given 
in terms of the Chebyshev polynomial of the first kind T/v 



P^{x) = 2A^T^{^), (3.10) 
where T]^{-^) = cos(A^cos~^(x/(2A)). The diagonal elements Xi of $ are also parameterized by 

Xj = 2Acos( ^' n). (3.11) 

Once the effective glueball superpotential is obtained, the free energy in the planar limit is then 
computed via 

^0 = ^1 W,f,iS)dS - ^SHogiA') (3.12) 

with the boundary condition that when there are other than quadratic terms in the tree level 
potential the free energy reduces to ()2.19p a.t gp = for all p 7^ 2. 

We will start with computing the exact effective superpotentials and the free energies for the 
quadratic, cubic and quartic tree level deformations. The results exactly agree with the matrix 
model calculations including the leading Veneziano-Yankielowicz term and all instanton corrections. 



3.1 Quadratic deformation 

First consider the quadratic deformation with n = 1 and g2 = m in ()3.5|) . In this case, using ()3.4p 
and (ITTTIl in gives 

W,s = mNA^ + —ul. (3.13) 
We then integrate in S and integrate out ui by minimizing 

W = mNA^ + ^uj - 2NS log(^) (3.14) 

with A and ui. This gives 

Wes = NS-NShgi^^), (3.15) 
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where Ai = [mls^y/^ is the scale of the M = 1 theory related to the scale A of the J\f = 2 theory 
by threshold matching of the gauge coupling running at energy scale m. ()3.15|) is exactly the 
Veneziano-Yankielowicz superpotential. Furthermore, putting in ()3.12|1 gives the same free 

energy ()2.19p as that obtained from the matrix model. In other words, the matrix integral in the 
planar limit of a zero dimensional bosonic matrix model is solved purely using integrating in and 
out and the Dijkgraaf-Vafa prescription. 



3.2 Cubic deformation 

Next let us consider the cubic potential with n = 2, g2 = m and = g in ()3.5|) . In this case, 
putting and ^T\} in (jSIHI) gives 

= {mN + 2guM' + —uj + (3.16) 



Minimizing 

with A and Ui gives 



,2 , 2 , 9 3 oAroT„_/^^ 



W = {mN + 2guM' + + - 2A^^log(-) (3.17) 



I^efr = NS- NS\og{^) + iV^log(l + ^u,) + ^ul + , (3.18) 

where 

mN m'^N^ B 



B = [-Ag'^N-'S + ^iV'^Wy -1 + 432^^ j (3.20) 

and Ai = (mA^)^/^. Expanding (j2.32j) and (j3.18p in g'^S/m^ gives exactly the same result including 
the leading Veneziano-Yankielowicz term and the instanton corrections to all order, 

IVefT = NS-NS\og{j,)-2Ns'—--NS{^f 

NS{^f - IQ2ANS{^Y - NSO{^f. (3.21) 



_^iVS(^)3 - 1024ArS(^)^ - NSO{^-^''' 
The free energy can also be computed putting ()3.18p or ()3.21|) in ()3.12|) . 



3.3 Quartic deformation 

Consider the quartic potential with g2 = rn, g4 = g and all other gp = in ()3.5|1 . Using ()3.4j) and 
(nmil in dSH), 

W., = mNA^ + ^ul + ^A^ + |n?A^ + (3.22) 
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Minimizing 
with A and ui gives 



W = mNA^ + + + + - log(^) (3.23) 



0(7 V ^ 

Putting (j3.24|) in (|3.23|) . we obtain the effective glueball superpotential 



w ATQ(^ ^TQ^ . Vl + 12(7g/m2-l , , ,^m^ A ^ 12(75 

^^'^ = ^^^2 - 12^) - "^^^^^^ 6AW^^ ) + ^12^ V ' + ^ ' ^'-''^ 

where Ai = (mA^)^/^. Making a series expansion of Wes given by p.39|l and ()3.25|1 again gives 
exactly the same result, 

= NS-NS\og{^) + lNS^-lNS{^^'' 



'Af^ 2 m2 2^ 



2 ' 



+^iV5(4)' - ?Ar5(4)' + iV5(9((4)')- (3.26) 
2 m"^ 4 m"^ m"^ 



Furthermore, putting (|3.25|) in (j3.12p and using the boundary condition discussed below (|3.12|) . the 
free energy matrix integral in the large N limit is computed using integrating in and out and the 
result indeed agrees with ()2.37|) . 



4 Conclusion 

We have shown that the Dijkgraaf-Vafa matrix model prescription with the proper normalization 
of the partition function and the scheme we presented gives the exact nonperturbative effective 
superpotential without a need for adding the Veneziano-Yankielowicz term by hand. Both the 
Dijkgraaf-Vafa matrix model prescription and the integrating in and out procedure give the same 
result for the full exact nonperturbative effective superpotential including the leading Veneziano- 
Yankielowicz term and all instanton corrections. It would be important to understand the physical 
relations between the two quite different structures, one based on the large N limit of zero di- 
mensional bosonic matrix models and the other based on the hypothesis of nonrenormalization 
and linearity principles. It is interesting to notice that because the parametrization of the fully 
degenerate Seiberg-Witten hyperelliptic curve was known in terms of the Chebyshev polynomials, 
the computation of the effective superpotentials on the integrating in side reduced to a simple sub- 
stitution of the parametrization of the diagonal elements of the adjoint scalar field into the gauge 
invariant coordinates and integrating in the glueball superfield. Consequently, one-cut matrix inte- 
grals are easily computed for any polynomial potential using integrating in and out. It would be 
interesting to find the polynomials that factorize the Seiberg-Witten curve for more general cases 
where the degeneration of the quantum moduli space is partial. Once these polynomials are found, 
doing multi-cut matrix integrals should not in principle be a difficult task using integrating in and 
out. Integrating in and out can thus be used as an alternative tool to do random matrix integrals. 
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